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We present optimal measuring strategies for the estima- 
tion of the entanglement of unknown two-qubit pure states 
and of the degree of mixing of unknown single-qubit mixed 
states, of which N identical copies are available. The most 
general measuring strategies are considered in both situations, 
to conclude in the first case that a local, although collective, 
measurement suffices to estimate entanglement, a non-local 
property, optimally. 

PACS Nos. 03.67.-a, 03.65.Bz 



I. INTRODUCTION 

Plenty of work has been performed in recent years on 
optimal quantum measurements, i.e. on measurements 
which provide the maximum possible information about 
an unknown quantum mechanical pure Jl]-^| or mixed || 
state, of which N identical copies are available. These 
works are focussed mainly on the determination of the 
unknown state as a whole, and consequently any of its 
properties is also estimated, although maybe not in an 
optimal way. 

On the other hand recent developments on the field 
of quantum information theory have stressed the impor- 
tance of the quantum correlations -or entanglement- dis- 
played by some states of composite systems. In the sim- 
plest of such composite systems, the two-qubit case, all 
non-local properties of pure states depend upon only one 
single parameter. Such non-local parameter is the only 
relevant quantity invariant under local unitary transfor- 
mations on each qubit and plays a central role in the 
quantification and optimal manipulation of entanglement 

In this work we analyze and solve the problem of opti- 
mally estimating the entanglement of an unknown pure 
state of two qubits. This problem has been independently 
addressed also by Sancho and Huelga in a recent work 
[fiJI, where only a restricted class of measuring strate- 
gies is considered. Here, on the contrary, we will con- 
sider most general quantum measurements on N iden- 
tical copies of the state. Their quality will be assessed 
through the gain of information they provide about the 
non-local parameter of the state. After presenting and 
proving the solution we will conclude that the optimal 
measuring strategies so defined are not equivalent to the 



ones used to fully reconstruct the unknown state. As a 
matter of fact, all information about some relative phase 
of the unknown state turns out to be irreversibly erased 
as the entanglement is estimated. 

Estimation of the degree of mixing of an unknown 
mixed state is a different but very much related topic 
that we shall also consider here. For the single-qubit case 
the amount of mixing is specified again by just one pa- 
rameter, the modulus of the corresponding Bloch vector, 
whereas in order to completely specifying the state two 
more parameters, namely the direction of the Bloch vec- 
tor, are also required. We shall show that in this case the 
optimal measuring strategy on any number N of qubits 
prepared in the same mixed state can be made compat- 
ible with the optimal estimation of the direction of its 
Bloch vector. 

Finally, we will show that a possible way of optimally 
determining the entanglement of an unknown, two-qubit 
pure state consists precisely in estimating, also optimally, 
the degree of mixture of any of its two reduced density 
matrices. Therefore, it turns out in this simple bipartite 
case that the optimal estimation of a non-local parameter 
can be done through a local measurement. 

The paper is structured as follows. Section [n] is de- 
voted to background material. We introduce a convenient 
parameterization of two-qubit pure states and consider 
their isotropic distribution. We also review some basic 
aspects on parameter estimation and on quantum mea- 



surements. In Section III we pose the problem of entan- 



glement estimation on firmer grounds and announce the 
main result of this paper: its optimal performance. Sec- 
tion [V, rather technical and that could well be skipped 
in a first reading, is devoted to the computation of some 
effective density matrix p( N '(b), an object which plays a 
central role in deriving the optimal strategy for estimat- 
ing entanglement. In Section [v] the N = 1,2,3 cases are 
presented in more detail in order to illustrate the general 
case. Optimal estimation of the degree of mixing is dis- 
cussed and solved in Section VI, and finally Section VII 
contains a discussion relating estimation of both entan- 
glement and mixing, and some concluding remarks. 



II. PRELIMINARIES 

We will consider here a two-party scenario. Alice and 
Bob will share the N copies of a completely unknown 
two-qubit pure state and their aim will be to obtain 
as much information as possible about its entanglement. 
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The sense in which the state is unknown, the mecha- 
nisms for extracting information from the system and the 
scheme for evaluating the extracted information will be 
briefly reviewed in what follows. 



A. Homogeneous distribution. 



case, for instance, the prior, unbiased density function for 
the parameter b is given by (||), so we have f(b) = 3b 2 . 
A generic measurement, allowing for the most general 
manipulation of the system, is represented by a resolution 
of the identity by means of a set of positive operators, 
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(4) 



All that is initially known about the state of each pair 
of qubits is that it is pure. This corresponds to the un- 
biased distribution on the Hilbert space H4 = H2 Ti-2 
of two qubits, that is, to the only probability distribu- 
tion invariant under arbitrary unitary transformations 
on . It is convenient to express the unknown state 
\ip) G H.2 ® Ti-2, which depends on six parameters, in its 
Schmidt-like decomposition 



IV>> = 
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b), (1) 



where the phase e %a , which is usually absorbed by one of 
the kets it goes with, has been left explicit. The non-local 
parameter b S [0, 1] characterizes the entanglement of 
Only for b = 1 is | ip) a product state \a) (g> \b), and 
thus unentangled. For b < 1 the state contains quantum 
correlations, b = corresponding to a maximally entan- 
gled state. Recall that this parameter is the modulus 
of the Bloch vector of the reduced density matrix pa on 
Alice's side, 



PA = tr B \tp){i>\ 



1 



\a)(d\ 



l-b, 



a)(- 



(2) 



and equivalently for ps- The other four parameters cor- 
respond to the two directions a and b of the Bloch vectors 
of pa and ps- Then, the unbiased distribution of pure 
states corresponds Q to the isotropic distribution of a 
in S 2 , b in S 2 , a in S 1 and the quadratic distribution of 
b in [0,1], 
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B. General measurements and information gain. 

The parties are thus provided with N copies of a pure 
state \tp) as in Eq. (Q), i.e. with the state lip)® 1 * 1 , and our 
aim is to construct the most informative measurement 
on the collective, 2A-qubit system for the estimation of 
the parameter b. The optimality criterion to be used 
is based on the Kullback or mutual information K[f, /] 
flit , a functional of two probability distributions /' and / 
that is interpreted as the gain of information in replacing 
the latter distribution with the former one |tL5fl. In our 



After the above positive operator valued measurement 
(POVM) has been performed, giving the outcome k with 
probability tr(M^ fc ^ p® N ), where p = we compute 

the posterior density function for b, f(b\k), through the 
Bayes formula 



f k (b) EE f(b\k) = 



P(k\b)f(b) 

P (k) ' 



where p(k) is given by 



p(k)= / dbf(b)p(k\b), 
Jo 



(5) 



(6) 



and the conditional probability of getting outcome k 
when the state's non-local parameter has value b, p(k\b), 
will be shown later. The gain of information resulting 
from obtaining the outcome k after the measurement is 
quantified by the Kullback information corresponding to 
the prior and posterior probability density functions 



K[f k ,f} = J dbf{b\k)\n( f -W-\ 



V f(b) J 



(7) 



This expression has to be averaged over all the possible 
outcomes of the measurement, so that the expected gain 
of information reads 

R[h,f}=^2p(k)K[f k ,f}, (8) 

k 

and using (|^) this expression can be written as 



K[f k J] = Y J J dbf(b)p(k\b) In 



P(k\b) \ 
p(k) J • 



(9) 



Let us notice here that the value of K[fk, /] in Eq. (Q) 
would remain unchanged if we decided to characterize 
the entanglement of \?p} by another parameter b = h(b) 
(where h(b) is any bijective function of the original pa- 
rameter b). Consequently, the gain of information we 
compute for b also applies to any of the measures of en- 
tanglement so far proposed, such as the entanglement of 
formation 



log 2 



log 2 



for the asymptotic regime, or the monotone [fiof 



1 - b 



(10) 



(11) 



for the single-copy case. 
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III. OPTIMAL MEASUREMENTS FOR 
ENTANGLEMENT ESTIMATION 

We are looking for a measurement of the form (Q) such 
that the expected gain of information (Q) is maximized. 
We will present and explain here and in Section ^ such 
optimal measurements, whereas their explicit construc- 



tion is mainly contained in Section IV 



A. Local and global strategies 

Before we proceed we comment on four classes of mea- 
surements Alice and Bob may consider in order to learn 
about b fjjj: 

• local measurements on only, say, Alice's side, i.e. on 
the N qubits supporting the local state N , would 
be the most restrictive class of the hierarchy; 

• uncorrelated bilocal -i.e. each party measuring on 
their local TV-qubit part independently- and 

• classically correlated bilocal -that is, with classical 
communication between Alice and Bob- measure- 
ments are two intermediate types of strategies; fi- 
nally, 

• global measurements on the 2N qubits constitute 
the most general case. 

Global measurements are in principle the most informa- 
tive ones. But as the parameter b which quantifies the 
entanglement of \ip), completely quantifies also the mix- 
ing of pa (and pb), it could well happen that local mea- 
surements, or bilocal on the two parties, optimal for the 
determination of the mixing, are as informative as the 
global ones with respect to entanglement. In fact, in re- 
ducing IV'XV'I to pa®Pb only the relative phase a is lost, 
the dependence on directions a and b and on the entan- 
glement b is preserved. We have found the optimal global 
and local measurement of b. The results obtained follow- 
ing the two strategies are the same, as we will discuss in 
Section VII , so all the extractable information about the 



entanglement is preserved under the partial trace oper- 
ation, and the four classes considered above turn out to 
be equivalent for entanglement estimation. 



B. Effective mixed state 

Notice that all the dependence on the measuring strat- 
egy (Q) in Eq. (||) is contained in the probability p(k\b) of 
outcome k conditioned to the entanglement of the state 
being some given b, 



P(k\b) 
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where the sum over the rest of parameters reflects the fact 
that we are only interested in the entanglement. This 
expression can also be written as 



p(k\b) = tr(M( fe V (JV) (^))> 
where the mixed state p( N \b) is 



P (N) (b) 



dd f db 
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(13) 



(14) 



Eq. (|l3| ) allows for an alternative interpretation to our 
problem: a 2iV-qubit mixed state p^ N \b) is drawn ran- 
domly with prior probability distribution f(b) = 36 2 and 
we want to determine it by estimating b. 

We will compute p(k\b) in the basis that diagonal- 
izes p( N \b), which will crucially turn out to be inde- 
pendent of b. Let us denote by Xi(b), ...,A m (6) the posi- 
tive eigenvalues of p( N \b), and with n%, ■-, n m their mul- 
tiplicity. From the normalization of ( |l4| ) the relation 
Sj=i n j^j = 1 follows. The sum n = J+ n j of multiplic- 
ities of (non- vanishing) eigenvalues equals the dimension 
of the space which supports \tp}(ip\® N . This is the sym- 
metric subspace of 7if N , and thus || 



(N + 2J)\ (N + 3)(N + 2)(N + 1) 
AH(2J)! 6 

With this notation Eq. (|l3|) reads 

rii ni+n 2 

p(k\b) = X 1 (b)J^M^ + X 2 (b) MgU 



(15) 



i—m-\-i 



,(*0 



(16) 



i— n — n m +1 



i=i 



By substituting this expression in (^|) and using the in- 
equality p6| 



[x\ + x-z) In 



x\ + x 2 

y\ + 2/2 



< X\ In ( — 

.Vi 
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,2/2 



(17) 



where Xi,yi > 0, along with the fact that the POVM is 
a resolution of the identity in the symmetric subspace of 



m 



i.e. 



j • 



it follows that the average gain 



of information is bounded by 



K[f k ,f}< 



dbf(b)J2^j(b)hx 



jdbf{b)\ 3 {b) 



(18) 
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C. Minimal most informative measuring strategy. 

The bound ( p"8|) can be minimally saturated through 
a measurement with m outcomes where each M^ k ' is 
the rife-dimensional projector over the subspace corre- 
sponding to the eigenvalue Aj, of p^ N \b), having then 
p(k\b) — nkXk(b). Therefore the construction of the op- 
timal measurement can be readily performed after the 
computation of the spectral decomposition of the state 
( |l4| ) , and this is done for an arbitrary N in the next Sec- 
tion. For a more detailed account of the N = 1,2,3 cases 
see Section where also the gain of information up to 
N = 80 has been computed explicitly. 

Notice also that there are other ways measuring strate- 
gies can be evaluated and, consequently, there is not a 
unique notion of optimality. For instance, in [Q-^) a guess 
for the unknown state is made depending on the outcome 
of the measurement, and then both guessed and unknown 
state are compared using the fidelity. It can be proved, 
following Ref. [H, that the optimal measurements pre- 
sented here, the most informative ones, are also optimal if 
we decide, alternatively, for a fidelity-like figure of merit 
satisfying some very general conditions |19[ . 



IV. COMPUTATION OF p <N) 

It has been shown that the spectrum of p( N > (6) deter- 
mines the maximal gain of information about b, whereas 
its eigenprojectors lead to the corresponding measuring 
strategy. Our next step will be the computation of the 
spectral decomposition of this effective mixed state. 

Let us rewrite the generic two-qubit pure state (Q) as 



U A ® U B {c.+\+) A < 
U A ®U B \^(b)), 



\+h 



(19) 



where c+ 
states 1+) 



c_ = \ ^r~, the single-qubit pure 



and 



constitute an or- 



thonormal basis in Alice's (Bob's) part -corresponding to 
some fixed direction in the Bloch sphere-, U a an d Ub are 
unitary transformations in each single-qubit space and 
\ip(b)) is a reference state. 

The state /)W (b) corresponds then to a Haar integral 
over the group SU(2) x SU(2), since it can be expressed 
as 



p {N \b)= dg{D(g)M(b)D(g)1) 



(20) 
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where the index g denotes the elements of the group G = 
SU(2) x SU(2), D(g) = U A ®U B is a \ x \ irreducible rep- 
resentation (irrep) of this group and M(b) — \ip(b)} (ip(b)\. 

A well-known result in group representation theory fol- 
lowing from Schur's lemma, the so-called orthogonality 



lemma, will be useful in the calculation of this integral. 
Consider a matrix A al3 (B) given by 



A^(B) = f dgD a (g)BD^(g), 

JgEG 



(21) 



where D a and D° are two unitary irreps of the group G. 
Then, 



Lemma 1 (orthogonality lemma): 

A al3 {B) = a{B)5 al3 I, 



(22) 



so A a/3 (B) is zero if the two representations are inequiv- 
alent and proportional to the identity if the two repre- 
sentations are equivalent. 

In order to benefit from this lemma we identify B with 
M(b)® N = \iP(b)}(iP(b)\® N and then consider the rele- 
vant irreps of SU(2) x SU{2) borne by the iV-fold tensor 



product of the 



-irrep of the group. These repre- 



sentations are the support of the state \ip(b)) 
next task is to recognize them. 

The state \ip(b))® N can be expanded as 



0A : 



and our 



\mf N = 4\++-++) A ®\-) l 



N-l 

-c I c_ 



N-2 2 
-C , C_ 



-c^- 3 c 3 



I A ^ I-/B 



-c N \— ) 



+ CiC 



N-l 



(23) 



where 



means that we have exactly the same vec- 



tor in the second subsystem. Notice that in the expres- 
sion above all the elements of the product basis {|iti)} 
of the local spaces Tlf N of Alice's and Bob's N qubits 
-i.e. K) = | + +... + +)> 2 ) = | + +... + -),---> 2 «) = 

| ... }- appear in the form \iii) A ® \ui) B . Notice, 

in addition, that if we denote by tot the sum of the 
third spin component of all spinors in each ket -i.e., 
for instance m T (\ + ++)) = 3/2, m T (|+H — )) = 1/2, 
tot (| — I — )) = —1/2, the terms multiplied by the 
same combination of the factors c+ and c_ have the same 
tot in A and B. The state ( ^3|) can thus also be expressed 
as 



®N _ N 



+c^- 1 c_ 



,.iV 



E 

_ N 

E 



\Ui) A <3 \Ui) B 



i;tnT=f-l 



\U„) , (X) \m 



(24) 
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We move now from the local spin basis {l^i)^} to the 
coupled one {l^i)^} in Alice's N qubits, and we also do 
the same in Bob's. The following lemma, that can be 
easily checked, will be useful here. 

Lemma 2: Let {|ej)} and be two orthonormal 

basis in C l , related by an orthogonal transformation O, 
so that \ ei ) = £\ Oijlfj), with O* = O, and O' 1 = 0*. 
Then, 



= £1/*)®!*: 



(25) 



i=i 



Now, notice that the unitary transformation relating 
the local basis and the coupled one is real (since all the 
Clebsch-Gordan coefficients are real) and that there is a 
conservation rule for the total third spin component (i.e. 
the Clebsch-Gordan coefficients that couple two states 
with third component rri\ and to.2 to a coupled state with 
third component to are proportional to 5 m ^ mi+m2 ). Then 
Eq. ( pi|) can be reexpressed, using the previous two facts 
and lemma 2, in the coupled basis as 



\i>(b)f N = 4 £ 

N 

c - E 

N 

2;mT — 



\ v i) A ® \ v i) B 



\Vi) A ®\Vi) B . 



(26) 



l\m T = — 4f 



(see the examples in next Section for more details). Wc 
note that the symmetry between the terms in A and in 
B allows us to derive (||^) from (^4|). 

Let us now have a closer look into Eq. (^). The 
term with coefficient corresponds simply to the state 
with a total spin j maximal in both Alice's and Bob's 
subsystem (i.e., j A = 3b 



y) and also maximal third 
spin component to, namely tua = tub — -y- We can 
thus write, with the notation \i A m A ) A <E> P b tob) b , = 
\v\) A ® \v\) B — \~? -f-) A ® I"^t)b - ^ n * s state belongs to a 
f <g> f -irrep of the group SU{2) x SU(2). The coefficient 
c+~ 1 c_ corresponds to all states with tra — m B — 



Apart from \v 2 ) = — 1)4 ® l~T"^)s' wm ch again 
belongs to the previous y ® y-irrep, the remaining N—l 



2^2 

kets, \v 3 ) ■ ■ ■ \v N+ x) have j A = Jb 



N 
2 



1, and thus 



belong to N-l different (but equivalent) (f—l)<8> (y-l)- 
irreps of the group. But since only the linear combination 
\v 3 ) + ■ • • + \vn + i) appears, the relevant irrep is just the 
symmetric combination of the latter N—l ones, which 
we will denote by {(y — 1) ® (t — ^)}sym, and which no 
longer decomposes as the product of two irreps of SU(2). 
The same applies for — 2) ® — 2) -irreps and so on. 



Thus, the space which supports the initial state can be 
decomposed in terms of irreps of SU(2) x SU(2) as 



N N 
~2~®~2~ 



N 



I] I y-1 

~ ' sym 

Nmod2 Nmod2 



(27) 



where Nmod 2 is equal to one for odd N and equal to 
zero for even N. It can be checked that this result agrees 
dimensionally with formula (|i~5|). 

The decomposition shown above in terms of the rel- 
evant irreps of the group SU(2) x SU{2) together with 
the orthogonality lemma can be used to solve the inte- 
gral in (^0|). As we have argued, when plugging (|2^) into 
( pp| ) the cross terms corresponding to inequivalent repre- 
sentations -such as |i>i)((u3| + .... + (vn+i\)- vanish as we 
integrate, while the terms within the same representation 
-such as lead to a contribution proportional to 

the identity in the subspace associated with the repre- 
sentation. So the state p^ N \b) is equal to 



p {N) {b) = \i{b)I^N + \2(b)I { (N_ 1 ^N_ 1 ) }sym - 

+ X m (b)I { Nmod 2 Nmod 2 y ^ 



(28) 



This is the spectral decomposition we are looking for, 
where {Xj} are the entanglement dependent eigenvalues 
of p( N \b), the trace of the identities giving the corre- 
sponding multiplicities {rij}. It is important to notice 
that, as it was mentioned before, the eigenspaces are in- 
dependent of b. 

The calculation of rijXj can now be readily performed 
from Eq. ( p6|) by computing the trace of the projection 
of \4>(b)) N into each relevant irrep. The determination 
of the spectrum of p( N ^(b) completes, as we have shown, 
the construction of the optimal measurement for the es- 
timation of the entanglement. In the next section some 
examples are studied in order to clarify the implementa- 
tion of the procedure. 



V. SOME EXAMPLES: THE N = 1,2,3 CASES 
AND BEYOND. 

In this section we will apply the procedure described 
above to obtain the optimal estimation of b when one, 
two and three identical copies of the initial state are at 
our disposal. 



A. TV 



The simplest case, N 
state written as in ( f[^ 
SU(2) x SU{2). From 
nality lemma as in (Eq) 



1, is now straightforward. The 
belongs to the | <£> \ irrep of 
we have, using the orthogo- 
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pW(b) = I dgD(g)M(b)D( g y = Ax (6)7. 



(29) 



The eigenvalue Ai(6) = \ is obtained by taking the trace 
in the expression above. The probability p(k\ b) (see (|l3|)) 
is independent of 6, so that p(k) = p(k\b) and the average 
Kullback information (^) vanishes. 

Consequently, no information whatsoever can be ob- 
tained about the entanglement of a completely unknown 
pure state if only one copy is at our disposal. 



B. N = 2 

For the N = 2 case the initial state has the form, from 

(ED ot dD» 

|VK6)} 82 = 4|++). 4 ®|.) B 

2 



+d|--) A ®|.) 



(30) 



Now, using lemma 2 and the conservation law mentioned 
before for the Clebsch-Gordan coefficients (cf. Eq. (p6[)), 
we can rewrite the state as 



|^(6)) 



02 



„2 II 



1) 



+ C+C _(| 1 0) A ®|.) B + | 0) A ®|.) S ) 



„2 II 



-1) 



I A w I-/B) (31) 

where for each party the coupled basis is related to the 
local one by means of an orthogonal transformation, as 
usual, 

l 1 l) = |++): I 1 



-1 = - - 



1 



|'0> = ^l + -> + |- 

l»o> = -L(| + -,-|-+) 



The state \il)(b))® 2 in |3l]) is supported then in the 1 <g> 1- 
and the <g> 0- irreps of SU(2) x SU(2), and now the 
application of lemma 1 gives for pW(b) 

p {2) (b) = Xi(b)I 19 i + Hb)I om . (33) 

We just need to pick up the contributions of (|3l]) to each 
irrep, that is the trace of the corresponding projections, 
to find that 

n 1 \ 1 (b) = (c% + <*<?_ + £) = 



n 2 A 2 (fe) = c 2 + c 2 



l-b 2 



(34) 



The optimal measurement (see Eq. (|l^)) then con- 
sists of two projectors onto the 1 <g> 1- and (g) 0-irreps of 
SU{2) <x> SU(2), with probabilities p(l\b) = niAi(6) = 
and p(2\b) = n 2 A 2 (6) = and from them 

p(l) = and p(2) = ^j. Finally the gain of information 
can be computed using (|) and it gives K — 0.0375 bits. 



C. N = 3 

The last case we want to discuss is N — 3. Starting 
now from (12(3) we have 



\mf=c\H) A ®\-) 



-cic_ 



+c+c^ 



+c°_ 



21 A 



(\*-h) A ®\-h 



+ * 



.3\ 
21 A 



IB- 



■IB 

(35) 



we observe that only contributions to the | <8> |- and to 
two different | <E> |- irreps of SU(2) x SU(2) appear. No- 
tice, in addition, that since in this expansion the contri- 
butions to ^ <g) ^ and to ^ (g> ^ only appear in a symmetric 
linear combination (i.e. |3 i)^ ® + |l i) ® |.) B and 

\')b + I* ~54® l-)s)' tne relevant irreps is pre- 
cisely a symmetric combination of the two latter ones, 
{ \ ® \ }sym- The orthogonality lemma gives now 



^ 3 )(6)=A 1 (6)/j 8 | + A 2 (6)/ { i 8 i }s 



(36) 



Finally, by collecting the traces of each projection of ( |35| ) 
onto each irreps we obtain 



mXxib) = (4 + c 4 +C 2 - + 44 + c 6 _) = 

n 2 A 2 (6) = 2 (c 4 + 4 + 44) = 



l + b 2 



(37) 



and thus the optimal measurement is composed by a 16- 
dimensional and a 4-dimensional projectors into the two 
irreps shown above, the corresponding probabilities being 
p{\\b) = ±±f- and p{2\b) = ±=£-. From them p{\) = § 
(32) andp(2) = ±, and the gain of information is of 0.084 bits. 



D. N > 3 

We have applied the same, general procedure to obtain 
the gain of information up to N = 80, as reported in Ta- 
ble I and Figure [l| We observe a logarithmic asymptotic 
dependence of the gain of information on the number N 
of available copies of which reads 

#«0.441og 2 iV (38) 

bits of information on b. 

VI. OPTIMAL ESTIMATION OF MIXING 

So far we have considered the most general measure- 
ment involving the whole space (H 2 ® of N copies 
of a two-qubit pure state. Now we are going to study 
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optimal local measurements for the estimation of its en- 
tanglement. Alice will perform a collective measurement 
over the N copies of the state pa in Eq. (0) at her dis- 
posal in order to estimate the parameter b. Consequently, 
we are also studying optimal strategies for estimating the 
degree of mixing of a singlc-qubit mixed state, when N 
copies are available. 

In order to study the latter with more generality we 
will consider a generic prior distribution f(b) for the de- 
gree of mixing while keeping an isotropic distribution in 
the Bloch vector direction a of the unknown mixed state, 
with 



that equivalent irreps, such as the N — 1 copies of the 
— 1)— irrep, obtain equal but independent contribu- 
tions. The space Tlf N , decomposed in terms of irreps of 
SU{2), is (see also Refs. § and @) 



N 

Y- 1 



N 



Nmod2 



Nmod2 



(44) 



The spectral decomposition of (b) is determined 
by application of the orthogonality lemma. Since equiv- 
alent irreps receive always the same contributions in the 
decomposition (J45| ), the corresponding eigenvalues are 
equal, so that fl4l|) reads 



s 2 



(39) 



A general measurement on the local composite system p^\b) = Af (6)/jv + A2 (b) ( 



I, 



supporting the state p^ consists of a resolution of 
the identity in the corresponding Hilbert space Hf N by 
means of positive operators . The gain of informa- 
tion is as in (^j), where now 

p(k\b)=tr(M^p A N \b)y (40) 

so that we need to compute the effective mixed state 



Pa 0) = 



gee 



dg (D{g)p A {b)D{ g y 



(41) 



where the integral is performed over the group G 
SU(2) and a single copy of the mixed state 



PA = U A p A {b)U\ 



(42) 



has been expressed, as before, in terms of a reference state 
PA(b) = (c+|+)(+| + c 2 _\— ){— I) and a unitary transfor- 
mation XJ a- The procedure to be followed is analogous 
to the previous one, the spectral decomposition of the 
state ( ffl| ) allowing us to build the optimal measurement. 

The density matrix pa{o)® n can be written -by using 
a straightforward modification of lemma 2 and the men- 
tioned properties of the Clebsh-Gordan coefficients- in 
terms of the coupled basis {l^i)^} as 



PA{b) 



>\A 



2!N-1) 2 



„2N 



— N 1 



i;rriT — — 



(43) 



Notice that the important role played before by the sym- 
metry between the kets in A and B (cf. Eq. (|26|)) is 
now played by the symmetry between the terms in the 
bra and in the ket. However we see that now there are 
no cross terms between inequivalent irreps of SU(2), and 



+ A^(6) [i Nmgd 2 + . . . + I n mod 2 ) . (45) 



This is, of course, simply what remains from Eq. ( |2q ) 
when Bob's subsystem is traced out, and we have in- 
cluded the whole derivation only for completeness. 

Eqs. (|l6- 18) still hold and therefore the optimal mea- 
surement for the degree of mixing b corresponds, for any 
isotropic distribution, to projections onto each of the sub- 
spaces associated with the eigenvalues {A^}. The gain of 
information is then given by the right hand side of Eq. 
(|l8|). Notice that both the number of outcomes and the 
corresponding probabilities p(k\b) — n^X^(b) are equal 
to the ones obtained before for entanglement estimation. 
In particular, it follows that there is no way to learn 
about the degree of mixture of an unknown mixed state 
if only one copy is available. 



VII. DISCUSSION AND CONCLUSIONS 

We have presented in this work an optimal strat- 
egy for the estimation of the entanglement of two-qubit 
pure states, when iV copies are available. Such opti- 
mal measurement is also minimal, in the sense that it 
consists of the minimum number of outcomes, namely 
A72 + 1 ((N + l)/2) outcomes for the even (odd) N- 
copy case. Most of the corresponding projectors are of 
dimension greater than one, and of course any further de- 
composition of them can be used in principle to obtain, 
simultaneously, some additional information about other 
properties of the unknown state, although our optimal 
POVM is not compatible with projecting onto states of 
the form Uj)i}® N as optimal POVM for state determina- 
tion do |2-p|, and they are thus less powerful for that 
purpose. 

An interesting particular case is when the initial state 
is a product one, i.e. 6=1. It can be seen that in this 
situation we have only the outcome corresponding to the 
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space of maximum spin, since niAi(l) = 1. Therefore 
if the outcome k, with k > 1, is obtained we can assure 
that the state is entangled. 

In the previous Section we have also been concerned 
with the optimal estimation of the degree of mixing. Our 
optimal measurement, again minimal, can be used, for in- 
stance, to quantify the degree of purity of states created 
by a preparation device whose polarization direction we 
ignore. Our strategy is actually complementary to the 
one aiming at revealing optimally the direction of polar- 
ization of the state jTj. As a matter of fact, the optimal 
POVM we have obtained is just a coarse graining of the 
one obtained in [0 for optimal estimation of mixed states, 
which turns out to reach also the optimal standards of 
direction estimation obtained in Q. Consequently, di- 
rection and modulus of the Bloch vector of an unknown 
mixed state can be optimally estimated simultaneously. 
Notice that this is not a frequent situation. If, instead, 
we would like to estimate the x, y and z components of 
the Bloch vector independently, we would have obtained 
incompatible optimal strategies (consider e.g. the N = 1 
case, where an optimal measurement for the component 
of the Bloch vector along direction h consists of a two 
outcome measurement projecting on that direction). 

Finally, we can argue that bilocal measurements, ei- 
ther uncorrelated or classically correlated, do not imply 
any improvement on the simpler, local ones for entangle- 
ment estimation. Once we get an outcome from Alice's 
local measurement we can compute Bob's effective state, 
and it is clear from Eq. fl2S| ) that his outcome will be the 
same as Alice's, so that no extra information on b will 
be obtained. We have also seen that the optimal global 
measurement on | - 0) 8JV is perfectly mimicked by a local 
one on p® N (or pf ), so that actua lly all four classes of 
measurements considered in Section [II A are equivalent. 
In fact, with hindsight, one can understand this result: 
local measurements are performed on the reduced den- 
sity matrix, which is obtained by a partial trace over the 
other subsystem. This operation erases the information 
contained in the parameters a and b of Eq. ([l]) . On the 
other hand the global measurement can be interpreted 
as being performed on the effective density matrix of Eq. 
(|l4|), where the same parameters have been integrated 
over. This operation erases the information contained in 
them too. 

It would be challenging to address the same question 
for bipartite mixed states, and for systems shared by 
more than two parties. Notice that in none of these cases 
optimal estimation of the non-local parameters would be 
possible by means of local (or even uncorrelated bilocal) 
measuring strategies. This is the case for mixed states 
because any given reduced density matrix pa may cor- 
respond to infinitely many mixed states p, with different 
degrees of entanglement, so that not even in the limit 
N — > oo can the entanglement of p be properly inferred 
from p®^. The mere existence of hidden non-local pa- 



rameters jl8j, that is of entanglement parameters that 
are erased during the partial trace operation, also pre- 
vents uncorrelated local strategies from being optimal for 
estimation of pure-state tripartite entanglement. 

To conclude, two-qubit pure-state entanglement, a 
quantum non-local property, can be optimally estimated 
by means of local, but collective, measurements. 
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N 


_ 

K 


1 





2 


0.03751 


3 


0.08397 


4 


0.13259 


5 


0.18059 


10 


0.39245 


20 


0.69639 


40 


1.07422 


60 


1.32005 


80 


1.50261 



TABLE I. Average gain of information K about b given TV 
copies of the state \rp). 



Average K 




FIG. 1. Average gain of information K about b given 
JV copies of the state \ip). The points represent the results 
obtained by the described optimal measurement, while the 
line shows the asymptotic behavior. 
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